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The temperature fluctuations of the Cosmic Microwave Background (CMB) are supposed to be 
distributed randomly in both magnitude and phase, foUowing to the simplest model of inflation. In 
this paper, we look at the odd and even multipoles of the spherical harmonic decomposition of the 
CMB, and the different characteristics of these, giving rise to a parity asymmetry. We compare the 
I ' even and odd multipoles in the CMB power spectrum, and also the even and odd mean angles. 

I , We find for the multipoles of the power spectrum, that there is power excess in odd multipoles, 

■ compared to even ones, meaning that we have a parity asymmetry. 

' Further, for the phases, we present a random walk for the mean angles, and find a significant 

, ^ , separation for even/odd mean angles, especially so for galactic coordinates. This is further tested 

^ ^1 and confirmed with a directional parity test, comparing the parity asymmetry in galactic and ecliptic 

coordinates. 
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INTRODUCTION 



Gaussianity and statistical homogeneity and isotropy of the CMB are essential for determining the most optimal 
values of the cosmological parameters in the framework of the ACDM concordance model. It is remarkable, that 
^ ' recently obtained data by WMAP, BOOMERANG, MAXIMA, CBI etc. clearly demonstrate a consistency between 
Qh, theoretically predicted and observed properties of the CMB anisotropy power spectrum([l[, Q and Q). However, 
Q ■ detailed investigations of the statistical pr^erties of the CMB signal reveal a significant deviation from Gaussianity 

" ■ ~ 1, i, i, a, iji, 0, [HI, [i3, [HJii, [H, in, [ii, [H, M), 



or it could even have primordial 



, and statistical isotropy of the CMB 
C/3 ' which can be related to the foreground residuals ([21|, [22]), systematic effects 

^1 origin (for instance the primordial magnetic field, texture or a non-trivial topology of the Universe). 

The important task of investigating the departure from statistical isotropy and Gaussianity of the CMB, requires 
, designed tests, which can clarify the possible sources of "contamination" in the results, derived from the CMB signal. 
Most of these tests are based on the CMB map, which require a good quality of data reduction, including a high quality 
of the calibrations, correction of the antenna beam and pixelization uncertainties and a reduction of the residuals of 
the point-like sources and diffusion foregrounds (and even the incorporation of the uncounted ones, like spinning dust 



\^ , or even the refracting debris of the Kuiper belt [23|). In contrast, recently proposed in [20|, [24|, a parity test only 



^y-^ ' requires the power spectrum, which is usually estimated with an accuracy significantly better than the CMB map. 

, Recently, ideas has been put forth, to explain the observed parity asymmetry, see for instance [2^ and [l^. Also 27] 
T-H ' has commented on the possibility of biased maps, which might also contribute to parity asymmetry. 

, The idea of our test is based on the analysis of the ratio g{l) of the powers stored in even and odd multipoles. For 
^ ■ a statistically homogeneous and isotropic random CMB, g{l) should be equally distributed around g{l) = 1. This 
assumption is applicable only, if the theoretical power spectrum C{1) of the primordial fluctuations, has a monotonic 
shape without any oscillations of the corresponding scale A/. The monotonic shape of the power spectrum is typical 
. for the simplest models of inflation, when the slow roll approximation is valid, until the epoch of reheating. However, 



in more complicated models (for instance in [28|), the slow roll regime of inflation could be broken in the vicinity of 
some points 4> C {4>i],i = 1,2., leading to spatial modulations of the power spectrum for adiabatic perturbations of 
the metric without departure from Gaussianity. In this case parity test can be used as indicator of non-trivial regimes 
of inflation. Additionally, we present a directional parity test G{1), similar to to the ratio g{l). This also tests the 
power stored in even and odd multipoles, but it is based on the directional statistic D{l,n) originally introduced by 
[29j. The directional parity test G{1) should also be equally distributed around 1. 

We would like to focus our attention on an extension of the parity test for a very broad range of multipole numbers 
I for the whole-sky CMB map, taking under consideration the phases (pi^m of the corresponding coefficients ai^m of the 
spherical harmonic decomposition. The phases (fn^rn carries much of the information about the angular distribution of 
the temperature anisotropy on the sky and possible correlations between different multipoles. Moreover, if the parity 
asymmetry of the power spectrum reflects possible departure from statistical isotropy of the CMB, the phases (j)i^m 



should have the same type of asymmetry. Previously, [30| and [3l| has investigated the phases, and has further found 
evidence for hemispherical asymmetry [32[ 



The outline of the paper is the following. In Section 1 we introduce the basic characteristics of the power and phases 
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asymmetry. For the phases we wih use the 'mean angle' approach, proposed in [17| for even and odd multipoles and 
show that there is strong departure from statistical isotropy. In section 2, we introduce the directional parity test. 
Section 3 takes a look on the alignment of mean angles, for low I. Section 4 is devoted to the random walk of the 
phases, in both galactic and ecliptic coordinates. In section 5 we present the results of the directional parity test, for 
ecliptic and galactic coordinates. Finally in section 6 we discuss the findings and draw our conclusions. 

ESTIMATORS OF THE PARITY ASYMMETRY OF AMPLITUDES AND PHASES. 

The temperature fluctuations on the sphere AT{8, (j>) can be decomposed into spherical harmonics in the following 
standard way: 

oo I 

AT(0,0)=^ ^ ai,™yi,™(0,0) (1) 

1=0 m=-l 

where n = {6, (j)) is a unit vector, pointing in the direction of the polar 9 and azimuthal 4> angles, ai^m is the coefficient 
of decomposition: ai^m = |ai,m| exp(i(/>i^m), with as the phase. 

Under the assumption of total gaussian randomness, as predicted by the simplest infiationary model, the am plit udes 
\ai.m\ follow the Rayleigh's Pdf and phases of ai^m are supposed to be evenly distributed in the range [0, 27r] [3J|. 
For any signals T{n) on the sphere, one can define the symmetric (T+(n)) and asymmetric (T~(n)) components 

r(n) = r+(n) + r-(n), (2) 

where 

, T(n)+r(-n) v-^ o,7r?, 
r+(n) - \ =E E a;,mC0s2(-)r,„(n), (3) 

I m— — l 

T^-(n) = "^^"^ '/^'"^ - E E «.™sin^(^)y,„(fi) (4) 

/ m— — l 

andr/„,(n) = (-iyy,„(-n). 

By definition, T+(n) and T~(n) are orthogonal, which means that the average over the whole sphere corresponds to 
(T+(n)r^(n)) = and that the power spectrum of T{h) is given by 

= ^E - C+il) + C-il) (5) 

m 

which is associated with power spectrum of even (C+(Z) — C{1) cos'^{^)) and odd [C^ [I) = C{1) fiiv?{^^)) multipoles 
respectively. For the concordance ACDM cosmological model with initial Gaussian adiabatic perturbations, we do not 
expect any features distinct between even and odd multipoles. However, taking the WMAP7 data as a platform for 
investigation of the parity of the CMB, there have been reported power contrast between even and odd multipoles of 
WMAP TT power spectrum (jsH], [ISl, 3^, [s^]) [S])- The corresponding estimator for "even and odd" asymmetry 

m = (6) 

The f;(/)-statistic is a cumulative statistic, but it is the most practical approach to testing for asymmetry between 
even and odd multipoles. The problem is, that we cannot make a comparison for a single 1, as either the odd or 
even component would be zero. An alternative could be to compare multipoles of a separation of 1 for instance, i.e. 
compare an odd multipole with its even neighbor, however we would have to decide what neighbor to compare to {I 
and 1 + 1, OT I and ^ — 1), and the choice would be an arbitrary one. In [s^ they performed a test for neighboring Cz's, 
comparing I and 1 — 1. We tested that method also, and found that it did not work as well as eq. [6] for our specific 
purposes, and thus we kept the (7(/)-statistics as introduced by [20^]. 

a At lowest multipoles (2 < / < 22), there is reported odd multipole preference (i.e. power excess in odd multipoles 



of the CMB power spectrum is given by 
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and deficit in even multipoles) ([35|, l20|, [36|, [37|). Note that due to rotational invariance of g(l), the corresponding 
parity asymmetry does not depend on the reference system of coordinates. However, additional information about the 
statistical properties of the temperature fluctuations is associated with the odd and even components of the phases 
(pi^rm which are rotationally non- invariant. This is why the analysis of the CMB phases (j)i^rn,, for different systems of 
coordinates seems to be extremely important for determining the origin of the CMB parity asymmetry. 
The phases carry information about the orientation of the spherical harmonic components, to make them fit the 
observed microwave radiation. In order to check the statistical anisotropy of the phases (j)i^m for / = odd and I ~ even, 
we will use the trigonometric moments: 

Si(0 - ^ E Ci(0 = ^ cos(</)/,„.) (7) 

The mean angle of the given multipole, I, is the arc+tangcnt to the ratio between the average sin and cos value for 
the multipole: 

e(0=arctan(^^) (8) 

The reason for the above procedures for finding the mean angle, is that we are dealing with circular data and thus 
a simple average over the phases is not enough. Instead we must convert the angles to points on the unit circle via 
sin /cos, and find the arithmetic mean of these points (Eq. [T]). This mean is then converted back to polar coordinates 
via arctan (Eq. [8]) to give us the mean angle for the considered multipole. Thus we have greatly compacted the 
information of the CMB phases, by creating a single mean angle, for each multipole. Further details about the 
procedures above can be found in [39|. 

We would expect, that the mean angles for the multipoles would distribute themselves evenly in the range from — tt to 
TT, as required by a Gaussian random field [40] . Moreover, for a statistical homogeneous and isotropic random Gaussian 
field generated by primordial perturbations, the uniformity of the distribution of phases leads to a uniformity of the 
distribution of mean angles 6^(0 for even and odd /. 

For testing the equivalence of the odd and even multipoles in the space of phases, we will use the following estimators: 
R^il) = j^{[rtil)f + [rtil)f}, r±G„,,)- XJsine±(0, r± = g cos e±(/), I > 3. (9) 

1=2 !=2 

where O"*" is the mean angle for even multipoles, G~ is the mean angle for the odd, and the corresponding sum is 
defined for n — 2k and n — 2k + 1, k ^ 1,2... for 0+ and O^. After simple algebra, Eq® gives us the following 
expression for the function R^: 



71,71 



l(k2k + [^s,^,k+i] + j^ E cos(e±(n)-e±K)); (lo) 
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From Eq(fTO|) one can see, that for non-correlated mean angles Q^jn), Q^{n') the last term can be estimated as /" 
for I » 2, and thus asymptotically R^{1) \ (see for details [41|). Actually, R"^ is well known in the theory of 
statistical analysis of circular data (see for instance, [i^l). The important generalization of this test, is that it will be 
applied for odd and even multipoles separately, describing a random walk of the mean angles in the space of phases. 
If the mean angles for given multipoles Q^{n),Q^{n') are close to each other (cos(0^(n) — 9^(n')) '-^ ±1), the 
second term in Eq pO|l is comparable with the first one. For the WMAP CMB map this effect is well known for the 
quadrupole and the octopole, as alignment. However, in the next section, we will show that the quadropole/octopole 
alignment is not unique. It is detected for AZ = 1 mean angles (even-odd multipole correlations) in the galactic system 
of coordinates, and especially, for AZ = 2 in ecliptic coordinates (see figd]). Thus the parity test of the mean angles 
G^(Z) needs to be taken under consideration as a complimentary test to the parity asymmetry of the power spectrum, 
in order to clarify the origin of these anomalies in the CMB. 
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THE COMPARISON OF THE DIRECTIONAL PARITY TESTS. 



Historically, investigation of the low multipole anomalies of the CMB signal by rotationally non-invariant (di- 
rectional) statistic was proposed in The estimator of the preferable direction n, which maximize the angular 
momentum(L) dispersion, can be defined as 

D{l,n) = {AT{n){{nLfAT{n)) - ^ m2|a,,„(n)|2 (11) 

m 

where the n is selected, so that the term (nL) has the maximum value. 

Actually, one can define the same estimator, as Ea ([TT]) . for the map of gradient in azimuthal direction AT^ = -g^AT{n) 
for a given system of coordinates with x-axis along the vector n. Then, based on the Eq ([TT|) . we can find the direction 
of maximal parity asymmetry of the CMB, for the considered directions, specified by n, just by replacing l{l -\- 1)C{1) 
with this new D{l,n) in Eq®, to get the following: 

In other words, an estimator precisely as Eq®, just using D{l,n) instead of C{1), for selected values of n. In this 
work, we only test for two n's, the one corresponding to galactic coordinates, and the one corresponding to ecliptic. 

For implementation of the parity test we use the Internal Linear Combination (ILC7) map. This is, according to 
the WMAP-team, the most precise data, when one is interested in m ultip oles / < 32, and it is indeed the very same 
data, that the WMAP-team uses for finding cosmological parameters [43| . 
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FIG. 1: Left panel. The parity parameters g{l), computed for the range of multipoles I < 40. The thick black line corresponds 
to WMAP7 power spectrum from the incomplete sky, the red line corresponds to the ILC7. Right panel. The same as left, but 
for 40 <l < 400. The red line represent the ILC7. 



TABLE I: Estimator for ILC 7 aligned muhipoles at 2 < Z < 40 

cos(e(0 - e(z-h 1)) cos(e(0-e(z-h2)) cos(e(0 - e(z -h i)) 

galactic coordinates ecliptic coordinates ecliptic coordinates 

1 = 2,3 0.9714 1 = 5,7 0.9986 / = 15, 16 0.9980 

= 18, 19 0.9947 I = 23,25 0.9995 I = 18, 19 0.9924 

= 28, 29 0.9978 / = 33, 35 0.9998 

= 33, 34 0.9477 / = 34, 36 0.9999 

= 36,37 0.9299 

= 38, 39 0.9485 



FIG. 2: Left panel. The mean angle for ILC7 phases in galactic coordinates, computed for the first 40 multipoles. Red dots 
indicates the odd multipoles, and blue dots are for even. A line between two points signify, that the two multipoles have a 
similar mean angle. Right Panel. The same as left, but for ecliptic coordinates. 




FIG. 3: Top panel. The parameter r~ (l) (the black dots) for odd multipoles and r^{l) (the red dots) for even ones in Galactic 
coordinates. Left graph is for the first 40 multipoles, right is for I up till 500. Bottom panels. The same as the top panels, but 
for ecliptic coordinates. 



STATISTIC OF THE PARITY PARAMETER AND THE MEAN ANGLES. 



For illustration of the method we select the WMAP 7 year data from http://lambda.gsfc.nasa.gov/, and following 



to the procedure above, we calculate the parity parameter g{l) for I < 400, and the mean angles for galactic, and then 
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FIG. 4: Left column. The parameters R~ (l) (top panel, the black line) for odd multipoles and R^{1) (top panel, the red line) 
for even ones versus / in galactic coordinates. Right column. The same as left, but for ecliptic coordinates. The black dash 
line is the asymptotic R'^{1 — even) = 0.5, the blue dash line is for _R~(Z = odd) = 2{i'~2) • Second and third from the top are 
for 1/R'{1) and l/i?+(0, at varying scales. 



for ecliptic coordinates, bearing in mind that the most statisticahy significant result for g{l) take place for 2 < I < 22. 
The results for the power spectrum are plotted in Fig [T] in comparison with parity parameter g[l) taken from the 
incomplete sky (the WMAP power spectrum). Since g{l) < 1 for all models at this range of multipoles, the power of 
the CMB C(2 < I < 40) has deficit of C{1 = even) modes in comparison with C{1 = odd) and vice versa. To clarify 
the origin of the parity asymmetry, we need to look closely at the complimentary test: the parity of mean angles. 

In the theory of circular statistics of phases it is well known, that for uniform and non-correlated (pi.m, the mean 
angles 6^(0 should be uniform and non-correlated as well. Similar to Eq([8l) we can define the mean angles £,^{1 C A) 
for a given range of multipoles, within some interval A. Namely 

^±^^A^ , -i /^ E;ca sin 9^(0 ^ 
\T,icA cose±(/) J 

The properties of the mean angle at any system of coordinates, depend on the correlation between the phases for a 
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given multipole I and the different m for this muhipole. In a separate paper, we wiU perform a general analysis of the 
most preferable direction of the phase correlations trough the mean angle statistics, trying to systematize the most 
peculiar directions. Here, for simplicity we will test the mean angle of the phases for even and odd in galactic 
and ecliptic coordinates. 

For the ILC 7 map in the Galactic coordinates at 2 < Z < 40 the odd multipoles (red dots in the figure [2]), have 
a mean angle of = 1.2238 rad in the range A = 2 — 20. The distribution of the even multipoles is characterized 
by the mean angle — 1.1919 rad. Note that without separation of the phases into odd and even one, we have 
for the mean angle } = 1.2080 rad if 2 < Z < 20. For the range 20 < ^ < 40 we have = 2.94286, C+ = 0.6783 
and ^ = 1.78893 rad. Thus, the cross-correlations between odd and even modes for this range A = 2 — 20 are 
given by K = cos(f+ - = 0.99949 and K = -0.639 for A = 20 - 40. For the ecliptic system of coordinates 
we have C~ = 1.06046, ^+ = 2.0799, ^ = 1-6152 rad. and K = 0.5238, where 2 < ^ < 20, and = -0.08576, 
^+ = -3.1331, C = -0.32940 rad. and K = -0.99556 for 20 < I < 40. Here we would like to direct attention to 
the remarkable similarity of the directions ^~ and for the ILC7 in Galactic coordinates, at the range 2 < I < 20. 
Based on uniformity of the distribution of the mean angles, the difference ^~ — should be uniformly distributed 
within the interval — 7r,7r. The probability that both these angles diflter by an angle A = — is given by 
P ~ A/(27r) ~ 0.0051. Note that this probability is close to the probability of getting the parity asymmetry of the 
power spectrum 0.0044 for the same range of multipoles, presented in [36J. Both of these directions are highly aligned 
with the mean angle of the first 19 multipoles (without separation into odd and even multipoles) in the Galactic 
system of coordinates. 

The alignment between ■f^, for 2 < / < 20 is not a unique feature for the mean angles of the ILC 7 phases. 
As it is seen from Fi^ there exists a coupling between even and odd mean angles, most obviously known for Z = 2 
and I = 3[ll|. In Table |T] we show the most interesting alignments between even and odd, odd-odd and even-even 
multipoles, for which the corresponding correlations exceeds the factor 90%. 

PEARSON'S RANDOM WALK IN THE DOMAIN OF MEAN ANGLES. 

In Fi^ we illustrate Pearson's random walk in Galactic and ecliptic coordinates, in the space of mean angles, in 
term of the parameters rf- versus r^, for the ILC 7 map at the range of multipoles 2 < / < 40. Note, that the random 
walk is of course dependent on the chosen coordinate system, so we expect different results from ecliptic and galactic 
coordinates. The upper limit I = 40 has been chosen, in order to avoid contribution from instrumental noise into the 
phases of the CMB signal. Also included, is a plot with / up till 500. The figure shows the separation of the even and 
odd mean angles in the space of phases. For ecliptic coordinates this tendency is quite weak, and the corresponding 
probability to get this effect from a random field is in the order of 1.5(7, where a = \/l/2 ~ 6.324. 

In galactic coordinates the statistical significance of the separation of even-odd modes exceeds ia threshold. More 
qualitatively the contribution from different multipoles into the even-odd random walk is illustrated in Fig|4l From 
this figure one can see, that a major contribution to the i?^(/)-parameter is given by the mean angles at 5 < Z < 30, 
where R~{1) exceeds the 2a threshold. The even multipoles, in contrast with the odd ones, are characterized by very 
small values of the i?+(^)-parameter, which indicate the presence of cross-correlations between even-even modes. 
It should be noted, that we expected values of R~{1) and R'^H) around 0.5 ( and thus 1/R ~ 2 for the i?^^-graph). 
That means that such very small values of R are as statistically unlikely, as very high ones. 

DIRECTIONAL G(/)-TEST. 

The result of the investigation of Pearson's random walk in the domain of the mean angles of the phases, clearly 
indicate significant differences in the separation between the even and odd multipoles, in the Galactic and ecliptic 
system of coordinates. The separation of the CMB phases for even and odd multipoles, is very clear established in the 
Galactic coordinates, and it is less significant in ecliptic one. To confirm these predictions from the Pearson's random 
walk statistic, we perform a parity test G{1) for the directional power D{1), given by Eq fTTI) . The results of the 
investigation are shown in Figl5]for the power D{1) and corresponding parity parameter G(/), for the case of galactic 
and ecliptic coordinates. As it follows from this figure, the dominance of odd modes over even ones occur in both system 
of coordinates, but the most significant departure from G 1 take place at Galactic coordinates. The comparison of 
G{1) test from Figl5]and g{l)-test from FiglT]show that for 2 < Z < 22 the significance of these two tests is practically 
the same. It would be worth it to note, that unlike the parity test g{l), defined on the power spectrum of the CMB 
signal, the directional G-test completely ignore the contribution of all m = modes, showing practically the same 
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FIG. 5: The left panel: the power _D(/)(in fiK^) for the galactic (black line) and ecliptic coordinates (red line). Right panel: 
the parity parameter G{1) for the galactic and ecliptic coordinates (black and red respectively). 



odd-even non-equality, as g{l). This feature is important for the determination of possible sources of contamination 
of the primordial CMB signal (which we believe, is parity neutral {g ^ 1)) by exotic foreground residuals, such as 
the Kuiper belt. Recently the simplest model of the dust emissivity from the Kuiper belt (KBOE) was discussed in 
[23} . where emission of the KBOE has been modeled as a pure quadrupole component with corresponding coefficients 
of the spherical harmonic decomposition a2m=o 7^ and a^^'^o — ^he ecliptic system of coordinates. In other 
words, the model only include contribution from the m — 0-modes. According to our G(Z)-test, this model for the 
foreground cannot contribute to the parity asymmetry, as the test exclude all m = 0- modes, and we still find a parity 
asymmetry. Thus, to fit the observational data, this model needs significant modification, which would be investigated 
in a separate paper. Generally speaking, implementation of directional statistics in the space of phases and the power 
D{1) seems to be very useful for determining the sources of parity asymmetry of the CMB. 



DISCUSSION 



In this paper we have presented tests of the parity asymmetry of the Cosmic Microwave Background. We have 
tested the power spectrum of the CMB for differences in power of the even and odd multipoles, and found a power 
excess in the odd multipoles, for the entire range considered. 

Additionally, we found peculiar alignments between mean angles for the odd and even multipoles. In galactic coordi- 
nates the separation between the Z-values, for significant alignments, was 1, and for ecliptic the separation was 2 for 
most of the aligned multipoles. 

Further, in order to investigate the parity asymmetry from the power spectrum, we looked at the corresponding mean 
angles, and took a closer look at the relationship between the even and odd multipoles in phase space via Pearson's 
Random Walk. We found, that there was a separation between the even and the odd multipoles, with strongest values 
for the galactic set of coordinates, where the statistical significance exceeded 3a. 

We also checked the i?(/)-values for various I, and found that the most significant contribution to the {I) -values 
happened for 5 < / < 30. The even multipoles are on the other hand characterized by small values of R^{1). 
Finally, we introduced the directional parity test, for various alignments of the coordinate system. We checked for 
galactic and ecliptic coordinates, and found, that the most significant departure from parity symmetry happened for 
the galactic system of coordinates, in accordance with the pearsons random walk test. 

In conclusion, it is obvious, that there exists a significant parity asymmetry between the even and the odd multipoles, 
especially so for galactic coordinates. A nongaussianity in the signal from WMAP can have both instrumental and 
cosmological origin, and the coming data from the Planck Survey will help us greatly in determining if we see 
instrumental error in WMAP or have hints of new, undiscovered physics. 



9 



ACKNOWLEDGMENTS 



We thank the anonymous referee for constructive comments and remarks. We acknowledge the use of the Le gacy 
Archive for Microwave Background Data Analysis (LAMBDA). Our data analysis made the use of HEALPix j45l |. 
This work is supported in part by Danmarks Grundforskningsfond, which allowed the establishment of the Danish 
Discovery Center. This work is supported by FNU grant 272-06-0417, 272-07-0528 and 21-04-0355. 



[1] G. Hinshaw, et al. Five- Year Wilkinson Microwave Anisotropy Probe (WMAP)Observations. The Astrophysical Journal 

Supplement, 180, 225-245, 2009. 
[2] N. Jarosik, et al. Three- Year Wilkinson Microwave Anisotropy Probe (WMAP) Observations. The Astrophystcal Journal 

Supplement, 170, 263-287, 2007. 
[3] C. L. Bennett, et al. First Year Wilkinson Microwave Anisotropy Probe (WMAP) Observations. The Astrophysical Journal 

Supplement, 148, 1-27, 2003. 

[4] L. -Y. Chiang, P. D. Naselsky, O. V. Verkhodanov and M. J. Way. Non-gaussianity of the derived maps from the first-year 

Wilkinson microwave anisotropy probe data. The Astrophysical Journal, 590, 65-68, 2003. 
[5] P. Coles, P. Dineen, J. Earl and D. Wright. Phase correlations in cosmic microwave background temperature maps. Mon. 

Not. R. Astron. Soc, 350, 989, 2004. 
[6] P. Vielva, E. Martnez-Gonzlez, R. B. Barreiro, J. L. Sanz and L. Cayn. Detection of Non-Gaussianity in the Wilkinson 

Microwave Anisotropy Probe First- Year Data Using Spherical Wavelets. The Astrophystcal Journal, 609, 22-34, 2004. 
[7] C.G. Park. Non-Gaussian signatures in the temperature fluctuations observed by the Wilkinson Microwave Anisotropy 

Probe. Mon. Not. R. Astron. Soc, 349, 313-320, 2004. 
[8] H. K. Eriksen, D. I. Novikov, P. B. Lilje, A. J. Banday and K. M. Gorski. Testing for Non-Gaussianity in the Wilkinson 

Microwave Anisotropy Probe Data: Minkowski Functionals and the Length of the Skeleton. The Astrophysical Journal, 

612, 64-80, 2004. 

[9] F. K. Hansen, P. Cabella, D. Marinucci and N. Vittorio. Asymmetries in the local curvature of the WMAP data. The 

Astrophysical Journal Letters, 607, L67-L70, 2004. 
[10] D. L. Larson and B. D. Wandelt. The hot and cold spots in the wilkinson microwave anisotropy probe data are not hot 

and cold enough. The Astrophysical Journal Letters, 613, L85-L88, 2004. 
[11] D.J. Schwarz, G.D. Starkman, D. Huterer and C.J. Copi. Is the Low-1 Microwave Background Cosmic?. Physical Review 

Letters, 93, 221301, 2004. 

[12] M. Cruz, M. Tucci, E. Martinez-Gonzales and P. Vielva. The non-Gaussian cold spot in Wilkinson Microwave Anisotropy 

Probe: significance, morphology and foreground contribution. Mon. Not. R. Astron. Soc, 369, 57-67, 2006. 
[13] A. de Oliveira-Costa, M. Tegmark, M. Zaldarriaga and Andrew Hamilton. Significance of the largest scale CMB fluctuations 

in WMAP. Physical Review D, 69, 063516, 2004. 
[14] C.J. Copi, D. Huterer, D.J. Schwarz and G.D. Starkman. Uncorrelated universe: Statistical anisotropy and the vanishing 

angular correlation function in WMAP years 13. Physical Review D, 75, 023507, 2007. 
[15] K. Land and J. Magueijo. Examination of Evidence for a Preferred Axis in the Cosmic Radiation Anisotropy. Physical 

Review Letters, 95, 071301, 2005. 
[16] A. Rakic and D.J. Schwarz. Correlating anomalies of the microwave sky. Physical Review D, 75, 103002, 2007. 

Mon. Not. R. Astron. Soc, 380, 466-478, 2007. raeth 
[17] P.D. Naselsky, L.Y. Chiang, P. Olesen and O. V. Verkhodanov. Primordial magnetic field and non-gaussianity of the 1-year 

Wilkinson microwave anisotropy probe (wmap) data. The Astrophysical Journal, 615, 45-54, 2004. 
[18] H. K. Eriksen, F. K. Hansen, A. J. Banday, K. M. Gorski and P. B. Lilje. Asymmetries in the Cosmic Microwave Background 

Anisotropy Field. The Astrophysical Journal, 605, 14-20, 2004. 
[19] J. Hoftuft, H.K. Eriksen, A.J. Banday, K.M. Gorski, F.K. Hansen and P.B. Lilje. Increasing evidence for hemispherical 

power asymmetry in the five-year WMAP data. arXiv:0903.1229v2, 2004. 
[20] J. Kim and P. Naselsky. Anomalous Parity Asymmetry of the Wilkinson Microwave Anisotropy Probe Power Spectrum 

Data at Low Multipoles. The Astrophysical Journal Letters, 714, L265-L267, 2010. 
[21] P. D. Naselsky, A. G. Doroshkevich and O. V. Verkhodanov. Phase cross-correlations of the wilkinson microwave anisotropy 

probe internal linear combination map and foregrounds. The Astrophysical Journal, 599, 53-56, 2003. 
[22] P. D. Naselsky, A. G. Doroshkevich and O. V. Verkhodanov. Cross-correlation of the phases of the CMB and foregrounds 

derived from the WMAP data. Mon. Not. R. Astron. Soc, 349, 695-704, 2004. 
[23] M. Maris, C. Burigana, A. Gruppuso, F. Finelli and J. M. Diego. Large Scale Traces of Solar System Cold Dust on CMB 

Anisotropics. arXiv: 1010.0830, 2010. 
[24] J. Kim and P. Naselsky. Anomalous parity asymmetry of WMAP power spectrum data at low multipoles: is it cosmological 

or systematics?. Physical Review D, 82, 2010. 
[25] A. R. Zhitnitsky P and CP Violation and New Thermalization Scenario in Heavy Ion Collisions. Nuclear Physics A, 853, 

2011. 

[26] F. R. Urban and A. R. Zhitnitsky The P-Parity Odd Universe, Dark Energy and QCD. \arXiv:101L24F5v l, 2010. 



10 



[27] C. J. Copi, D. Huterer, D. J. Schwarz and G.D. Starkman Bias in low-multipole CMB reconstructions I arXiv: 1 1 03. 3505\ 
2011. 

[28] D. Polarski and A. A. Starobinsky. Spectra of perturbations produced by double inflation with an intermediate matter- 
dominated stage. Nuclear Physics B, 385, 623-650, 1992. 

[29] M. Tegmark and A. de Oliveira-Costa High resolution foreground cleaned CMB map from WMAP. Physical Review D, 
68, 123523, 2003 

[30] C. Raeth, G. E. MorfiU, G. Rossmanith, A. J. Banday and K. M. Gorski. Model-Independent Test for Scale-Dependent 
Non-Gaussianities in the Cosmic Microwave Background. Physical Review Letters, 102, 131301, 2009. 

[31] C. Raeth, A. J. Banday, G. Rossmanith, H. Modest, R. Suetterlin, K. M. Gorski, J. Delabrouille and G. E. MorfiU. Scale- 
dependent non-Gaussianities in the WMAP data as identified by using surrogates and scaling indices. I arXiv:101 2. 2985\ 
2010. 

[32] G. Rossmanith, C. Raeth, A. J. Banday and G MorflU. Non-Gaussian signatures in the flve-year WMAP data as identifled 

with isotropic scaling indices. Mon. Not. R. Astron. Soc, 399, 1921-1933, 2009. 
[33] C. Raeth, P. Schuecker and A. J. Banday. A scaling index analysis of the Wilkinson Microwave Anisotropy Probe three-year 

data: signatures of non-Gaussianities and asymmetries in the cosmic microwave background 
[34] J. M. Bardeen, J. R. Bond, N. Kaiser and A. S. Szalay. The Statistics of Peaks of Gaussian Random Fields. The 

Astrophysical Journal, 304, 15-61, 1986. 
[35] K. Land and J. Magueijo. Is the Universe odd?. Physical Review D, 72, 101302, 2005. 

[36] J. Kim and P. Naselsky. Anomalous parity asymmetry of WMAP power spectrum data at low multpoles: is it cosmological 

or systematics?. Physical Review D, 82, 063002, 2010. 
[37] A. Gruppuso, et. al. New constraints on Parity Symmetry from a re-analysis of the WMAP-7 low resolution power spectra. 

\arXiv:1006.T979^ 2, 2010. 

[38] C. L. Bennett and et al. Seven- Year Wilkinson Microwave Anisotropy Probe (WMAP) Observations: Power Spectra and 

WMAP-Derived Parameters. arXiv:1001.4758vl, 2010. 
[39] N. I. Fisher Statistical analysis of Circular Data Cambridge University Press, UK, 1993. 

[40] P. G. Ferreira, J. Magueijo and K. M. Gorski. Evidence for non-gaussianity in the cobe dmr 4 year sky maps. The 

Astrophysical Journal Letters, 503, L1-L4, 1998. 
[41] P. Naselsky, L. -Y. Chiang, P. Olesen and I. Novikov. Statistics of phase correlations as a test for non-Gaussianity of the 

CMB maps. Physical review D, 72, 063512, 2005. 
[42] K.V.Mardia and P.E.Jump Directional statistics Wiley, 2000. 

[43] D. Larson, et al. Seven- Year Wilkinson Microwave Anisotropy Probe (WMAP) Observations: Power Spectra and WMAP- 
Derived Parameters. arXiv: 1001.4635u2, 2010. 

[44] D. A. Varshalovich, A. N. Moskalev and V. K. Khersonskii. Quantum Theory of Angular Momentum. World Scientific, 
1988. 

[45] K. M. Gorski, E. Hivon, A. J. Banday, B. D. Wandelt, F. K. Hansen, M. Reinecke and M. Bartelmann. HEALPix: A 
Framework for High-Resolution Discretization and Fast Analysis of Data Distributed on the Sphere . The Astrophysical 
Journal, 622, 759-771, 2005 



